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We study surface and bulk properties of porous films produced by a model in which particles incide 
perpendicularly to a substrate, interact with deposited neighbors in its trajectory, and aggregate 
laterally with probability of order a at each position. The model generalizes ballistic-like models by 
allowing attachment to particles below the outer surface. For small values of a, a crossover from 
uncorrelated deposition (UD) to correlated growth is observed. Simulations are performed in 1 -|- 1 
and 2-1-1 dimensions. Extrapolation of effective exponents and comparison of roughness distributions 
confirm Kardar-Parisi-Zhang roughening of the outer surface for a > 0. A scaling approach for 
small a predicts crossover times as and local height fluctuations as at the crossover, 

independently of substrate dimension. These relations are different from all previously studied 
models with crossovers from UD to correlated growth due to subsurface aggregation, which reduces 
scaling exponents. The same approach predicts the porosity and average pore height scaling as 
and , respectively, in good agreement with simulation results in 1 -|- 1 and 2-1-1 dimensions. 

These results may be useful to modeling samples with desired porosity and long pores. 

PACS numbers: 81.15.Aa, 05.40.-a, 68.35.Ct, 68.55.-a 


I. INTRODUCTION 

A widely studied model for growth of porous films 
is ballistic deposition (BD)b in which the particles in¬ 
cide perperdicularly to the substrate and aggregate at 
the first contact with the deposit^i^. BD was originally 
proposed to describe sedimentary rock formationii and 
was extended to model thin film growth and related 
systems by considering other aggregation mechanisms, 
non-colimated particle flux, or polydispersivity of par¬ 
ticle sizem^—. Most works on the ballistic-like models 
address the scaling features of the outer surface of the 
deposits, particularly for the connections with Kardar- 
Parisi-Zhang (KPZ)i^ roughening. Some works also con¬ 
nect the surface growth dynamics and the bulk properties 
of the porous depositsidii^li. This is an essential step 
for proposing models of porous materials, which have a 
large variety of technological applications, frequently in 
the form of thin film o^^i^° . 

Some ballistic-like models are in a class of competitive 
growth models in which uncorrelated deposition (UD) is 
obtained for a certain value of a parameter and, near that 
value, a crossovers in kinetic roughening is observed^d^. 
In the simplest model, particle aggregation follows the 
BD (UD) rule with probability p {1 — p). It was already 
studied numerically^!^ and with scaling approaches^^Si. 
For small p, there is an enhancement of characteris¬ 
tic times of the correlated kinetics (p = 1) by a fac¬ 
tor p~^ and of the outer surface roughness by a factor 
for a recent discussion of this topic, see Ref.^. 
These features extend to other ballistic-like models with 
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crossovers to UDi^ and are related to the lateral aggrega¬ 
tion. In case of surface relaxation after aggregation, the 
exponents in those relations are larger, corresponding to 
longer crossovers for small 

A renewed interest in these competitive models was 
recently observed, with a focus on the limitations of scal¬ 
ing relations or with an emphasis on the properties of 
porous media. RefSf^^S^ discussed the deviations from 
the dominant scaling of surface roughness at low p, which 
is essential for a quantitative characterization of surface 
properties in a variety of growth conditions. The effect of 
relaxation after collision of incident and aggregated par¬ 
ticles was considered in Ref.>”, also with a focus on surface 
properties. Refsj ^^i^° considered the effect of a stickness 
parameter on the aggregation of the incident particles, 
which may attach to neighboring particles located be¬ 
low the outer surface. Simulations in 1 -1- 1 dimensions 
produced deposits with porosity ranging from very small 
values to approximately 70% and suggested non-KPZ be¬ 
havior in one of the models^S. 

The first aim of this paper is to study surface and bulk 
properties of the model proposed in Ref^ combining a 
systematic analysis of simulation data and a scaling ap¬ 
proach for small values of the stickness parameter. From 
the extrapolation of saturation roughness and relaxation 
times, we show that the model has KPZ exponents in 
1 -b 1 dimensions. Comparison of roughness distributions 
confirms KPZ scaling in 2-1-1 dimensions, thus ruling out 
the proposal of non-KPZ exponents. In the limit of small 
stickness parameter a, the scaling approach shows that 
the crossover time and the roughness scale as and 

respectively, for all substrate dimensions. These 
results show a shortened crossover when compared to all 
previously studied models with an UD componeni^ i^^i^'^ , 
which is a consequence of subsurface aggregation. The 
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same approach predicts porosity and pore height scaling 
as and respectively. These predictions will 

be confirmed numerically in 1 + 1 and 2 + 1 dimensions. 
The approach can be extended to the model introduced 
in Refi^, with the same crossover exponents due to the 
similar subsurface aggregation conditions. 

The rest of this work is organized as follows. In Sec. |lTl 
we present the sticky particle deposition model. In Sec. 
cni we analyze the surface roughness scaling of simulated 
deposits in 1 + 1 dimensions. In Sec. we present a 
scaling approach that relates surface properties to the 
stickness parameter, and confirm those predictions with 
numerical simulations. In Sec. El we extend the scal¬ 
ing approach to relate the porosity and the average pore 
height with the stickness parameter, again with support 
from numerical simulations. In Sec. ED we show that 
the main results extend to 2 +1 dimensions. In Sec. EID 
we summarize our results and conclusions. 


II. STICKY PARTICLE DEPOSITION MODEL 

In all models discussed in this paper, square particles 
of size c are sequentially released on randomly chosen 
columns of a one-dimensional discretized substrate of lat¬ 
eral size Lc and move in a direction perpendicular to the 
substrate. Here, L is the number of columns. The time 
interval for deposition of one layer of atoms (L atoms) is 
T. Thus, at time t, the number of deposited layers is t/r. 

The model proposed in Refj^ is hereafter called sticky 
particle deposition (SPD). In each site of the trajectory 
of the incident particle, it interacts with particles in near¬ 
est neighbor (NN) sites at the same layer (same height 
above the deposit) and particles in next neareast neigh¬ 
bor (NNN) sites at the layer immediately below it. This 
interaction is represented by a probabilistic rule of aggre¬ 
gation at its current position. 

The probability of aggregation to each neighbor is 

Pag = 7 TTri ) ( 1 ) 

(r/c) 

where a is the stickness parameter, r is the distance 
between the centers of the particles (r/c = 1 for NN, 
rlc = \f2 for NNN), and n is an exponent related to the 
nature of the interaction. In Refi^, the cases n = 2 and 
n = 6 were respectively called Coulomb-type and van 
der Walls-type interactions, with most results being pre¬ 
sented for the former. Here we will restrict the analysis 
to the case n = 2, in which aggregation to NN and NNN 
have probabilities a and a/2, respectively. 

Fig. [T]helps to understand the aggregation rules of the 
SPD model and the differences from other ballistic-like 
models. We first recall the rules of BD and of the next 
nearest neighbor BD (BDNNN) model^i^i^. 

In BD, aggregation occurs at the first contact with a 
NN occupied site: particle A at position 2, particle B at 
position 6 in Fig. [TJ In BDNNN, aggregation occurs at 
the first contact with a NNN occupied site: position I 




FIG. 1: Illustration of the rules of particle aggregation in the 
SPD model. Deposited particles are gray and black sqnares, 
the latter indicating those particles in the outer surface. In¬ 
cident particles are indicated as squares marked A and B and 
circles numbered from 1 to 7 indicate possible aggregation 
positions of those particles. 


for particle A and position 5 for particle B in Fig. [TJ In 
both cases, the incident particle interacts only with the 
particles at the top of each column, which are highlighted 
in III The set of top particles is called the outer surface 
of the deposit. 

In the SPD model, particle A may aggregate at lattice 
sites marked with circles labeled 1 to 4, and particle B 
may aggregate at lattice sites marked with circles labeled 
5 to 7. 

First, consider the trajectory of particle A. In position 
I, two aggregation trials are executed due to the inter¬ 
action with two NNN occupied sites; the probability of 
aggregation in each trial is a/2 . If it does not aggregate 
there, it moves to the position labeled 2, in which three 
aggregation trials are executed: two for interactions with 
the NN at the same height (probability a for each one) 
and one for interaction with the NNN in the layer be¬ 
low, at the left (probability a/2). If the particle does 
not aggregate at position 2, then it moves to position 3 
and may aggregate there with probability a due to the 
interaction with the NN at the left. If aggregation does 
not occur in position 3, the incident particle will perma¬ 
nently aggregate at position 4, which is the top of the 
incidence column. Relaxation to neighboring columns is 
not allowed. 

Now we consider the trajectory of particle B. In po¬ 
sition 5, two aggregation trials are executed, each one 
with probability a/2 (due to interactions with two occu¬ 
pied NNN sites). If the particle does not aggregate there, 
it moves to position 6, in which three aggregation trials 
are executed: two for the interactions with the lateral NN 
(probability a for each trial) and one for the interaction 
with the NNN in the layer below (probability a/2). If no 
aggregation trial is successful at position 6, the particle 





















3 


moves to position 7 and aggregates there. 

In contrast to other ballistic-like models (e. g. BD and 
BDNNN), the SPD model allows subsurface aggregation. 
In Fig. [T1 position 3 is an example of subsurface position: 
it is not allowed in BD, nor in BDNNN, nor in any model 
of solid-on-solid deposition (which prescribe aggregation 
at the top of each column). 

In all cases, note that the interaction of an incident 
particle with an aggregated one is possible in two steps: 
the first one when they are NNN (larger distance), the 
second one when they are NN (smaller distance). It rep¬ 
resents two possibilities of aggregation in the ingoing part 
of the trajectory of the incident particle. If the aggrega¬ 
tion trials are not accepted, then the incident particle 
moves to a lower position. In this situation, this particle 
is in an outgoing trajectory respectively to those aggre¬ 
gated particles. For this reason, no aggregation trial is 
executed with a NNN aggregated particle in the layer 
above the current position of the incident particle. For 
instance, when particle A is at position 3 (third layer of 
the deposit), we do not execute aggregation trials with 
the black NNN sites at the fourth layer in Fig. [TJ 

The SPD model resembles the model introduced in 
Ref.— and the slippery BD model (SBD) proposed in 
Ref—, both studied in three-dimensional deposits (the 
latter with line seeds perpendicular to a flat inactive sur¬ 
face). Most of our simulation work is in 1 -|-1 dimensions, 
similarly to Ref—, but in Sec. IVIl we show that the main 
results are also valid in three-dimensional samples. 

For simplicity, in the following sections we consider 
unit values of the lattice constant and of the time of 
deposition of a layer: c = 1, r = 1. 


III. KINETIC ROUGHENING 


The outer surface roughness is defined as 


W{L, t) = ({h — K) 


-T-\ 2 


1/2 


( 2 ) 


where h is the height of the top particle of each column, 
the overbars indicate spatial averages, and the angular 
brackets indicate configurational averages. In systems 
with normal (in opposition to anomalous) scaling, the 
roughness follows Family-Vicsek (FV) scaling^! as 


W{L,t)^L^f(^^y (3) 

where a is the roughness exponent, is a relaxation 
time, and / is a scaling function. In long times {t ^ tx), 
f const, so that W saturates as 

Wsat - AL^, (4) 

where A is a model-dependent constant. The saturation 
time tx scales as 




BL^ 


( 5 ) 



a 2.5 3 3.5 4 


log,o(L) 


FIG. 2: (Color online) (a) Time evolution of the surface 
roughness in the SPD model with a = 0.5 (red solid curve), 
a = 0.1 (green dashed curve), and 0.025 (blue dotted curve). 
The dashed line has slope 1/3 of KPZ scaling, (b) Saturation 
roughness as a function of the lattice size for a = 0.5 (red 
squares) and 0.1 (green triangles). The dashed line has slope 
1/2 of KPZ scaling. 


where 0 is the dynamic exponent and B is another model- 
dependent constant. The roughness for < -C tx scales as 

W (6) 

with 13 = ajz and another model-dependent constant C. 

Fig. [5^ shows the surface roughness evolution of the 
SPD model for three values of o in L = 1024. 

For short times, there is a crossover from an initial 
regime of rapid roughness increase to a second regime in 
which it increases slower. For small a, the first regime 
is mainly of UD and the slope of the log W x log t plot 
is near 1/2. For a > 0.1, lateral aggregation is frequent, 
thus the roughness at short times is larger than that of 
UD (e. g. a = 0.5 in Fig. [2^). It is difficult to find a 
pure UD regime in this case and to estimate the crossover 
time with accuracy. 

After this transient, the growth regime begins, with 
apparent power law scaling of W [Eq. ([S])]. It is difficult 
to distinguish the different curves for small a in Fig. [^; 
this will be explained by the scaling approach of Sec. |IV1 
The slope of those curves are near 1/3, suggesting KPZ 
scaling. 

At long times, there is an increase in the saturation 
roughness as a decreases. 

Fig. shows the saturation roughness as a function 
of lattice size L for two values of a. They seem to be con¬ 
sistent with the KPZ exponent a = 0.5. However, linear 
fits of those plots give slopes slightly smaller than 0.5, 
similarly to what was found in Ref—. For this reason, a 
systematic extrapolation of those results is necessary to 
decide whether the roughness scaling is KPZ or not. We 
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FIG. 3: (Color online) Effective roughness exponents as a 
function of L~^ for a = 0.1 (red triangles) with A = 0.72 
and a = 0.5 (blue squares) with A = 0.52. 

proceed by using the same methods of in 

which roughness scaling of various ballistic-like models 
was studied. 

Effective roughness exponents are defined as 

=-k2-^ <"> 

Assuming that the saturation roughness has scaling cor¬ 
rections as Wsat ^ L°'{ao + where oq and oi 

are constants, we expect fv a + 02 ^“^, where 02 is 
another constant. 

Fig. [3] shows effective exponents as a function of L~^ 
for a = 0.5 and a = 0.1, respectively using A = 0.52 
and A = 0.72. These values of A provide the best linear 
fits of the ttL X data for each stickness parameter. 
The asymptotic (L — 00 ) estimates from those fits are 
a = 0.496 ± 0.015 and a = 0.51 ± 0.02, respectively. 

We estimate the dynamical exponent z using the 
method proposed in Refi^. For each lattice size L, a 
characteristic time to is defined as 

Wito) = kW,at. (8) 

with fc < 1. The FV relation ([3]) shows that to is propor¬ 
tional to tx for fixed k, thus to ^ L^. Effective dynamical 
exponents are defined as 

_ ln[to (T)Ao (A/2)] 

« = -ta2-■ 

Figs. 0^ and HId show zl for a = 0.5 and a = 0.1, 
respectively, obtained with fc = 0.8. In both cases, the 
exponents oscillate near 0 = 1.5, suggesting that finite- 
size corrections are very small. 

The estimates of a and z are in very good agreement 
with KPZ exponents a = 1/2 and z = 3/2, which is a 


FIG. 4: (Color online) Effective dynamical exponents as a 
function of l/L for (a) a — 0.1 and (b) a = 0.5. 

strong numerical evidence that this model is in the KPZ 
class in 1 -b 1 dimensions. 

A universal scaling is expected in the SPD model be¬ 
cause there is no change in its symmetries as the stick¬ 
ness parameter changes. In other words, the correspond¬ 
ing hydrodynamic growth equation may have coefficients 
dependent on the parameter a, but the leading spatial 
derivatives will be the same^ii^i. Due to the lateral ag¬ 
gregation and consequent excess growth velocity, KPZ 
scaling is expected for any a > 0. 

Previous works on ballistic-like models^ii^^i^zii have 
already shown that systematic extrapolation of finite-size 
or finite-time data are necessary to avoid crossover ef¬ 
fects. As highlighted in Ref4i, this is a consequence 
of the large fluctuations in height increments, typical of 
those models. 

Crossovers and finite-size corrections probably are the 
reasons for the deviations from KPZ scaling observed 
in Ref/^. This may also be inferred by comparison 
with finite-size BD data from Refi^. Refj^ suggests 
a ~ 0.42 for o = 1, while Ref.— gives effective expo¬ 
nents 0.40 < aL < 0.45 for BD in the same range of L. 
Moreover, the growth exponents /3 in lattice sizes from 
L = 256 to L = 1024 for a = 1, shown in Ref—, are 
very near the corresponding estimates for BD in Ref— 
(considering minimum linear correlation coefficient 0.999 
in the growth region). 


IV. SCALING FOR SMALL STICKNESS 
PARAMETER 

A. Scaling approach 

For small a, lateral aggregation is unprobable, thus 
most particles aggregate at the top of the column of in- 
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FIG. 5: Two neighboring columns after some time of random 
growth. Deposited particles are gray and black squares, the 
latter indicating particles at the outer surface. The incident 
particle is the empty square and circles indicate possible ag¬ 
gregation positions of this particle. 

cidence. At short times, the roughness increases as^ 

( 10 ) 

After a crossover time tc, KPZ scaling appears. Our first 
step is to relate tc to a, for a <C 1 . 

A typical configuration of two neighboring columns in 
UD is illustrated in Fig. [S] It has a height difference 

( 11 ) 

because their heights increase without correlations. If 5h 
is large, then a new particle inciding at the right column 
in Fig. [5] may aggregate at a number of positions of 
order 5h, as indicated by the circles. This means that 
the number of aggregation trials is of order Sh and the 
aggregation probability for each trial is of order a. Thus, 
the probability of no lateral aggregation after those trials 
(i. e. aggregation at the top of the column) is Ptop 
(1 — a)^^. The probability that some lateral aggregation 
occurs is, consequently. 

Plat = 1 - Ptop 1 - (1 - a6h. ( 12 ) 

The latter approximation requires Piat aSh 1, which 
will be confirmed below. 

The average time for a lateral aggregation event at a 
given column is Uat ^ t/Plat- Lateral aggregation imme¬ 
diately creates correlations between the heights of neigh¬ 
boring columns, thus the crossover time is tc ~ Uat- Eqs. 
OT and (ini) at the crossover {5h = She, t = tc) give 

/ 1 \l/2 

Shetc^^"^ ^ ^ {aShc) (13) 

Thus, height fluctuations at the crossover scale as 

She ~ (14) 


and the crossover time scales as 

tc ^ (15) 

These results confirm that Piat ^ aSh < aShc <C 1, as 
the approximation in Eq. (HU requires. 

The amplitudes of the saturation roughness [Eq. ([3|)] 
and of the relaxation time [Eq. ([S])] scale as She and tc, 
respectively. Following the exponent convention intro¬ 
duced by Horowitz and Albano^i^, we have 


A - a-^ 

(5 = 1/3, 

(16) 

B^a-y 

y = 2/3. 

(17) 


These results are valid in any spatial dimension because 
UD properties are not dimension-dependent. 

The scaling exponents in Eqs. (HU) and (II3 differ from 
those obtained in other competitive models with ballistic- 
like aggregation with probability p and UD with proba¬ 
bility 1—p; in those systems, S = 1/2 and y = 1— In 
solid-on-solid models with crossovers from UD to corre¬ 
lated growth, the exponents are also different: (5 = 1 and 
y = 2. The shorter crossover of the SPD model is due to 
subsurface aggregation, which provides a large number 
of opportunities (of order She) for lateral aggregation of 
the incident particle (Fig. [5]). On the other hand, the 
relation S = y/2 obtained in other competitive models is 
also obeyed here because it is solely related to UD scaling 
(see e. g. the discussion in Ref;^). 

B. Numerical test 

We performed simulations of the SPD in L = 1024 and 
small values of a, from 0.1 to 0.0125, until the steady 
states (roughness saturation). The saturation roughness 
Wsat and the characteristic times to were calculated fol¬ 
lowing the same lines of Sec. IIIII 

Figs, and|n )3 show to and Wsat, respectively, as a 
function of the stickness parameter a. Since they were 
measured for constant L, they are expected to scale as 
the amplitudes B [Eq. (HI)] and A [Eq. ([J])], respectively. 
Fits of the data for a < 0.1 give exponents y ~ 0.70 [Eq. 
(fT71) and Fig. H^] and S ~ 0.27 [Eq. (flBl) and Fig. [Bjo]. 

The estimate of y is in good agreement with the the¬ 
oretical prediction of Eq. m- However, the estimate 
of S shows a discrepancy of ~ 20% from the theoreti¬ 
cal prediction of Eq. (HU). Note that the fits in Figs. 
[ 6 ^ and [ 6)3 considered 0.05 < a > 0.0125, which are not 
very small values of a, thus deviations are expected, par¬ 
ticularly in the smaller exponent (S). Unfortunately, it 
is very difficult to obtain accurate estimates for smaller 
values of a because relaxation times become very large 
and roughness fluctuations also increase. Using smaller 
system sizes is also inappropriate because it enhances 
crossover effects. 
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FIG. 6: (Color online) (a) Characteristic time to and (b) sat¬ 
uration roughness Waat in size L = 1024 as a function of 
the stickness parameter. Solid lines are fits of the data for 
a < 0.05. 


Refi^ estimated the crossover times for L = 512 and 
obtained tc ^ which is significanly different 

from the theoretical prediction in Eq. (HSJ. However, 
measuring reliable crossover times is a difficult task, as 
explained in Sec. m On the other hand, the same work 
shows that the saturation time for L = 1024 scales as 
g-o.7±o.03^ which is in good agreement with our estimate. 

The scaling of the amplitude C in Eq. (0) can be 
predicted along the same lines of Refs.— for other 
competitive models: 

, i = 5-yp. (18) 

Eor the SPD model in 1 -|- 1 dimensions, we obtain 7 = 
1 /9. This very small exponent gives a very slow variation 
of C with the stickness parameter. It explains the small 
distance between the curves for different values of a in 
Fig. Hi. 


C. Scaling in a related model 

In Refj^, a model similar to the SPD was introduced. 
The particles incide vertically and, at each site of its tra¬ 
jectory with a NN occupied site, it may aggregate with 
probability p. Otherwise, the particle moves down one 
site. If no lateral aggregation occurs, the particle aggre¬ 
gates at the top of the column of incidence. 

In Fig. H particle A may aggregate to positions labeled 
2, 3, and 4. In positions 2 and 3, aggregation trials have 
probability p. If the particle does not aggregate at one of 
those points, it moves to position 4 and aggregates there. 
Particle B may aggregate at position 6 with probability 
p, otherwise it moves to position 7 and aggregates there. 


FIG. 7: (Color online) Regions of size 48 x 48 (in lattice units) 
of samples grown with stickness parameters (a) 0.1, (b) 0.01, 
(c) 0.001, and (d) 0.0001. 


For small p, most lateral aggregation trials are rejected, 
thus UD dominates. Large local height fluctuations ap¬ 
pear, similarly to Fig. [5l The increase of the local height 
difference Sh and the probability of lateral aggregation 
Plat are given by Eqs. m and (HU, with a replaced by 
p. Thus, the same reasoning of Sec. IIV Al leads to the 
same scaling relations of the SPD model with a replaced 
by p. 

In the notation of Refj^, exponents a' = 1/3 and 
z' = 2/3 are predicted by our scaling approach. The nu¬ 
merical estimates of that work, a' ~ 0.25 and z' ~ 0.77, 
differ from those predictions, probably because they were 
obtained by data collapse methods that do not account 
for scaling corrections. 


V. POROSITY AND PORE HEIGHT SCALING 

For a = 1, the samples have large porosity $ Ri 70%. 
When a decreases, $ decreases because UD creates no 
holes. Fig. H shows regions of some samples obtained 
with small stickness parameters. The porosity decrease is 
accompanied by the formation of longer pores extended 
in the vertical direction. This is a consequence of the 
increase of the height fluctuation She before a lateral ag¬ 
gregation event [Eq. (fHl) ] . 

The number of deposited layers necessary to attain a 
steady state value of $ is relatively small, typically of the 
order of tc in Eq. (Ha. This is expected because pores 
are narrow, even in pure BD, thus porosity depends only 
on short wavelength height fluctuations, which saturate 
at short times (in the absence of scaling anomaly^). 

Our scaling approach can be used to predict the de¬ 
pendence of the porosity 4> and the average pore height 
on the parameter a, as follows. 
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During the time interval tc between two lateral aggre¬ 
gation events, the number of particles deposited at the 
top of a given column is approximately (note that we 
are still using unit lattice constant and unit deposition 
time of a layer, c = 1 and r = 1). The size of a long pore 
produced by the lateral aggregation is 5hc [Eq. (fT4)l ]. 
Consequently, for small a, the porosity (pore volume di¬ 
vided by total volume) is expected to scale as 


$ 


She 

tc She 


1 / 3 ^ 


(19) 


This is valid in the limit of very small o, in which She 

tc- 

The small exponent in Eq. m explains why a large 
decrease of a leads only to mild reduction of porosity. 
This is remarkably illustrated in Eig. [3 in which a varies 
three orders of magnitude, while the porosity decreases 
from $ Ri 0.39 to 0.044, i. e. changes by a factor smaller 
than 10. 

The porosity scaling in the SPD also differs from other 
competitive models involving ballistic-type aggregation. 
Examples are the bidisperse ballistic depositions^ and the 
BD-UD competitive model, in which $ ^ (p is the 
probability of the ballistic-like component). 

We simulated the SPD in size L = 1024 for small values 
of a in order to measure the porosity between times tj = 
5000 and tp = 10000. In all cases, tj is much larger 
than the crossover time and tp is much smaller than the 
relaxation time tx ■ 

Fig. shows the porosity as a function of the stick- 
ness parameter. The linear fit for 10“"^ < a < 10“^ gives 
$ ^ in excellent agreement with Eq. (jT^ . Al¬ 

though these values of a are very small, the corresponding 
values of and of 4) are not very small. Thus, scaling 
corrections are particularly weak in this case. 

For very small a, the pores are long and isolated, as il¬ 
lustrated in Fig. [3i- The average pore height is expected 
to scale as Eq. da, because a pore is formed only when 
a lateral aggregation event occurs. However, for $ ^ 0.1 
or larger, many pores occupy two or more neighboring 
columns. This can be observed in Figs. [IK, Eb, and [It. 

Here we define pore height as the vertical distance be¬ 
tween the aggregation position and the top of the inci¬ 
dence column in any lateral aggregation event. Its av¬ 
erage value, {hp), is taken over all lateral aggregation 
events between tj and in 10^ different samples. For 
small $, pores are isolated, thus (hp) is a reliable ap¬ 
proximation of the average pore height, and is expected 
to scale as Eq. m- For $ not too small, some pores 
occupy two or more neighboring columns, and all these 
columns contribute to {hp) (each one had a lateral ag¬ 
gregation event). 

Fig. |Bb shows (hp) as a function of a. The slope of 
that log-log plot evolves from —0.15 for 10“^ < a < 10“^ 
to —0.28 for 10“^ < a < 10“"^. The latter is 16% smaller 
than the theoretically predicted value —1/3 [Eq. (jH])], 
which indicates the presence of large scaling corrections. 



logio(a) 


FIG. 8: (Color online) (a) Porosity as a function of the stick- 
ness parameter. The solid line is a linear fit of the data for 
a < 0.01. (b) Average pore height as a function of the stick- 
ness parameter. Dashed lines have slopes —0.15 (right) and 
-0.28 (left). 


Ref.— measured the porosity of samples with 0.1 < a < 
1, with results in qualitative agreement with ours. How¬ 
ever, the low porosity scaling was not addressed there. 

Ref— suggests that the porosity scales with p (equiv¬ 
alent to a) and with the lattice size L. The latter is ex¬ 
pected only as vanishing corrections, since porosity does 
not depend on long wavelength fluctuations. This ex¬ 
plains the small (effective) exponents a and c obtained 
in that work. On the other hand, Ref— estimates the 
long-time scaling on p with exponent b = 0.22, which is 
to be compared with the theoretical prediction 1/3. The 
discrepancy is probably related to the use of data collapse 
methods. 


VI. SPD MODEL IN 2 + 1 DIMENSIONS 

The aim of this section is to show that the main fea¬ 
tures of the SPD model in 1 -|- 1 dimensions can be ex¬ 
tended to 2 -|- 1 dimensions, namely the KPZ roughening 
of the outer surface and the porosity scaling derived by 
the superuniversal approach of Sec. lYl 

The aggregation rules of the SPD model have to be 
extended in this case. First, NN interactions are consid¬ 
ered in two substrate directions, with a total of four NN 
in the same height. Secondly, NNN interactions appear 
with aggregated particles in the same height (four neigh¬ 
bors) and with particles at the level immediately below 
(four neighbors). 

Roughness scaling of ballistic-like models usually show 
large correction a^^'^^i^^ . An alternative to search for the 
universality class of a given model is the comparison of 
scaled roughness distributions of relatively small systems 







FIG. 9: (Color online) Scaled roughness distributions in the 
steady states of the SPD model with a = Q.l (squares) and 
of the RSOS model (solid curve) in 2 + 1 dimensions, with 
L = 256. 

because the finite-size corrections of those quantities are 
much smaller—!^. 

We simulated the SPD model with a = 0.1 in sub¬ 
strates of lateral size L — 256 up to the steady state 
(roughness saturation). In this regime, the square rough¬ 
ness W 2 = h'^ — h of several configurations is measured. 
P {w 2 ) is the probability density of the square roughness 
of a given configuration to lie in the range [w 2 , W 2 + dw 2 ] ■ 
This quantity is expected to scale as 

= , 20 ) 

where a = \J (^2^) — {w 2 )'^ is the rms fluctuation of W 2 
and 'I' is a universal function^^— . 

Fig. [9] shows the scaled roughness distribution of the 
SPD model and the distribution of the restricted solid- 
on-solid (RSOS) model^ in substrate size L = 256. The 
latter is a well known representative of the KPZ class 
and its roughness distributions have negligible finite-size 
effects^. The excellent collapse of the curves in Fig. 0 
is striking evidence that the SPD model also belongs to 
the KPZ class in 2 -I- 1 dimensions. 

We also simulated the SPD model in size L = 1024 
for small values of a and measured the porosity between 
times ti = 1000 and tp = 2000. We observe that the 
porosity is larger than in the 1 + 1-dimensional samples 
for the same value of a. For instance, for a = 0.1, the 
porosity exceeds 50%. This is a consequence of the larger 
number of interactions of the incident particle with NN 
and NNN in 2 -|- 1 dimensions, which facilitates lateral 
aggregation. 

Fig. (TU] shows the porosity as a function of a, for low 
values of that parameter. The linear fit for 10“^ < a < 


FIG. 10: (Color online) Porosity as a function of the stickness 
parameter of the SPD model in 2 -f 1 dimensions. The solid 
line is a linear fit of the data for a < 0.01. 


10 ^ gives $ ~ a which is also in excellent agree¬ 
ment with Eq. (IISl). This supports the extension of the 
scaling approach of Sec. |TV]to 2-1-1 dimensions. 

An important consequence of this scaling approach is 
to facilitate the design of samples with the desired values 
of porosity and elongate pores. However, one has to take 
care with the fluctuations in the value of $ in the first 
layers of the deposit, typically produced at t < tc- 


VII. CONCLUSION 

We studied surface and bulk properties of porous de¬ 
posits produced by a model proposed in Refj^, in sub¬ 
strates with one and two dimensions. The model shows 
a crossover from uncorrelated to correlated growth for 
small values of the stickness parameter a. 

In 1-1-1 dimensions, a systematic analysis of simulation 
data for saturation roughness and relaxation times shows 
that the model belongs to the KPZ class. Finite-size 
corrections explain the previous claim of deviations from 
KPZ scaling. In 2 -f 1 dimensions, KPZ roughening is 
confirmed by comparison of roughness distributions. 

A scaling approach for small values of a is proposed 
to relate the crossover time and the local height fluctua¬ 
tions with that parameter, respectively giving exponents 
—2/3 and —1/3. These results are consequence of the 
UD properties, thus they do not depend on the spatial 
dimension. Numerical results confirm these predictions. 
The crossover exponents are smaller than those of other 
competitive models that consider aggregation only at the 
outer surface^iS^iiH. 

The same approach predicts the porosity scaling as 
a}/^, which is in good agreement with simulation results 
in 1-1-1 and 2 + 1 dimensions. This result is important for 
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using the model to produce porous samples representa¬ 
tive of real materials. This may also help to model sam¬ 
ples with desired porosity and pore height, particularly 
for the possibility of controlling the scaling properties by 
changing the kinetics of subsurface aggregation. 
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